

  *********************************************
  *******		CSES - COUNTRY FE		*******
  *********************************************
  
clear all	
cd "${data}"	
use "data1.dta", clear 


gen observationlimit_CSES = . // in order to limit the bivariate model to the identical sample with the full model.
replace observationlimit_CSES = 1 if daysafter!=. & IMD5012_!=. & election_!=. & country_!=. & respondent!=. 

* bivariate model 1
eststo CSES_countryFE0: mixed IMD3007_ IMD5012_ i.country_ if observationlimit_CSES==1 || election_: || respondent:

*base model
eststo CSES_countryFE1: mixed IMD3007_ daysafter IMD5012_ i.country_ if observationlimit_CSES==1 || election_: || respondent: 


eststo CSES_countryFE2: mixed c.IMD3007_ c.daysafter##c.IMD5012_ i.country_ || election_: || respondent: IMD5012_  



cd "${tables}"	

  esttab CSES_countryFE0 CSES_countryFE1 CSES_countryFE2 using tableL1.tex , b(3) se(3) nogap transform(ln*: exp(2*@) 2*exp(2*@)) label mlabels ("Model 1" "Model 2" "Model 3")
  